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Abstract 

I show that the discrepancy on the critical flavor number of fermions Nf for 
the appearance of the infrared fixed point based on the t'Hooft anomaly 
matching condition and derived from the Schrodinger functional method 
(Nj ~ 9) and the experimental analysis of the JLab group using Bjorken 
sum rule and GDH sum rule, and our lattice simulation (Nj ~ 3) could be 
resolved by assuming the topological structure of the infrared fixed point is 
not that of U(l) 3 but that of G2 with triality automorphism of octonions 
which appear in the product of quaternions. 

The agreement of the infrared fixed point of the running coupling mea- 
sured in lattice simulations with use of the quaternion real condition and the 
prediction of the BLM renormalization theory might be due to the boundary 
condition on S 3 x R manifold of quaternion. The form factor of a proton 
measured via Ward identity through the difference of the inverse propagator 
at momentum p+q/2 and at p-q/2 agrees with the phenomenological dipole 
fit. 

Keywords: Infrared fixed point, Critical flavor number, Chiral symmetry 
breaking, Quaternion, Triality, G2 symmetry 



1. Introduction 

The effective number of flavors of fermions in QCD is important in two 
aspects. The first is the anomaly cancellation, i.e. since axial vector cur- 
rent Ward identity is not automatic, it is necessary to seek an anomaly 
cancellation mechanism to assure the renormalizability of the field theory. 
The second is from the hadron phenomenology based on the perturbative 
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QCD(pQCD). Based on the quenched lattice simulation data, Banks and 
Zaks[6| presented a conjecture, that when the energy becomes low and the 
QCD effective coupling becomes large, the beta function crosses a zero and 
an infrared fixed point appears when infrared singularity is suppressed by 
large number of quark pair creations. The flavour number above which a 
fixed point appears is called critical flavor number N c . The search of N c in 



lattice simulation was performed in Schrodinger functional method [17| and 



Appelquist et al[14| claimed that it is around 10. 

Although running coupling below QCD scale Aqcd is difficult to interpret 
in pQCD, the running coupling at high energy can in principle extended to 
low energy via effective charge method^, §], and how to interplet the critical 
flavor number obtained in Schrodinger functional method is important. 

One can calculate effective coupling in lattice and in Dyson-Schwinger 
equation(DSE), via ghost-gluon coupling and via quark-gluon coupling. In 
the lattice simulation of the ghost-gluon coupling, we observed strong fluctu- 
ation of the ghost propagator and slight suppression of the infrared effective 
coupling. In the DSE approach of quenched QCD, the effective ghost-gluon 
vertex in Landau gauge is calculated as, 

a{p ) = a(fi )■ 



where D \p 2 ) = = ^ (p 2 + m 2 ). 

The term m 2 suppresses the fluctuation and the 'scaling solution', i.e. 
assume that the gluon propagator z(p 2 ) ~ (p 2 )~ KA and the ghost propagator 
G(p 2 ) ~ (p 2 )~ KC with k = Kq — — yields the running coupling which is 
not infrared suppressed (3]. Although the qualitative behavior of the momen 



turn dependence of the effective coupling of DSE and lattice simulation |26 
agree, I think the suppression of fluctuation via mass the term is artificial, 
and that the fluctuation of the ghost propagator is one of causes of the sup- 
pression of the infrared divergence of the ghost-gluon coupling. The DSE 
approach suffers from momentum cut-off, which could circumvented through 
the Pinch Technique [81]. 

The effective quark-gluon coupling in Coulomb gauge showed clear signal 
of infrared fixed point, and the result is compatible with the experimental 
values extracted from the polarized electron-proton scattering data 



via Bjorken sum rule[34[ and Drell-Hearn-Gerasimov sum rule |35l. |36| . The 




data of JLab group suggest infrared fixed point of a s (0) ~ 3. From their 
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fitted function of the effective coupling, the beta function obtained from the 
fitted line does not cross zero. The similar beta function was obtained also 
in the ADS/QCD approach Q. 

In the infrared region, color confinement is a characteristic feature and 
we measured Kugo-Oiima parameter 22] in quenched [23( and unquenched 



gauge configur at ions |2 51. |27jj. We observed that the Kugo-Ojima criterion 
is consistent with unquenched lattice simulation but about 20% deviation 
occurs in quenched simulation. This qualitative difference suggests strong 
correlation between gauge field and the quark field. Violation of the BRST 
symmetry in quenched calculation is observed also in DSE approach [lo|, [Til, 



12, 0]. In 27 1, I showed reduction of fluctuation in effective coupling and 
effective mass of quarks by choosing a gauge that the quarks on domain 
walls are correlated by self-duaj_gauge field (instanton). 



In the book of E. Cartan [45], a vector field like gauge field is produced 
from spinor fields as Pliicker coodinates and the spinor is expressed in quater- 
nion basis. I analyze in this paper, the lattice results of the Coulomb gauge 
quark propagator using quaternion basis and study consequences of the topo- 
logical structure of the manifold. A product of quaternion makes an octonion 
and the octonion has the symmetry of exceptional Lie group G2. It has the 
symmetry, which is called triality. It has the effect of triplicate flavor num- 
bers. Although in real world, quarks are massive and the system is not 

conformal, I expect proximity of our system to the conformal window. 

he 

In the infrared, when the Compton wave length of a quark where 

mc 2 

m is the effective mass of a quark is large (~ 0.3GeV), the proximity of 
other quarks in a baryon will change the boundary condition of quarks. I 
study the charge form factor of a proton using quarks in quaternion basis 
and modifying boundary condition of quarks. It makes a correction like final 
stated interaction and make the result gauge dependent. 

This work is organized as follows. In sect. 2, I summarize lattice data 
related to the IR fluctuation, and in sect. 3, I summarize the Schrodinger 
functional method. In sect. 4 anomaly matching condition is reviewed and in 
sect. 5 the octonion and anomaly matching condition are explained. I explain 
our method of gauge fixing of the gluon and domain wall fermion system 
in sect. 6. The form factor of the proton calculated by this gauge fixing is 
presented in sect. 7. 

I discuss whether our finding of infrared fixed point in Nf = 2 + 1 lat- 
tice simulation is incompatible with the estimates of Schrodinger functional 
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method and give conclusions in sect. 8. 



2. Unquenched lattice simulations and QCD effective coupling 

We performed lattice simulation of the effective coupling in Landau gauge 
from ghost-gluon coupling and found that it is close to the decoupling solu- 
tion of the DSE. We found that the fluctuation of the ghost propagator in the 
infrared is large and the deviation of the Kugo-Ojima parameter in quenched 
simulation from the value expected from the BRST symmetry is due prob- 
ably to this fluctuation. While in unquenched simulation, the fluctuation is 
suppressed and Kugo-Ojima parameter turned out to be consistent with the 
theory, and I guess the fluctuation is an essential property of the IR QCD. 

In the case of effective coupling, we calculated the product of the ghost 
propagator and the gluon propagator and observed infrared suppression due 
to suppression of the ghost propagator. The ghost propagator in the Landau 
gauge has color antisymmetric part whose expectation value is but its norm 
becomes comparable to that of color diagonal part in the infrared (28[. This 
kind of problem is absent in the Coulomb gauge. The lattice data calculated 
by using the gauge configuration of MILC |29|] and 16 3 x 32 x 16 domain wall 
fermion of RBC/UKQCD [30| suggest an infrared fixed point of a s (0) ~ 3 in 
the case of Coulomb gauge [271] . 



The formula of Brodsky-Lepage-Mackenzie(BLM) scheme [371. |38| strongly 
suggests a presence of infrared fixed point of a s (0) = ir and the effective 
coupling has walking like behavior near infrared. In the BLM scheme, ob- 
servables A and B are related in terms of the effective charge as 

®a(Qa) = azB{QB)(l + r A /B 1 ) (1) 

7T 

and the coefficient va/b is so chosen that the result does not depend on the 
number of flavors Nf. 

I show in Fig. (TJ the effective coupling of JLab experiment, and that of 
the lattice simulation of the quark gluon coupling. 

In our lattice simulation of Nf = 2 + 1, a s (q) ~ 1.7 at q ~ 0.6GeV and 
the deviation of the running coupling from 2-loop perturbative calculation is 



significant below q ~ 3GeV due to A 2 condensates [39|, |40(. I think the non- 
perturbative effect is significant in the infrared and the Schrodinger functional 
method do not disprove the infrared fixed point at Nf = 2 + 1. In 19|, |20|, the 



origin of appearance of the infrared fixed point is attributed to the acquiring 
of mass of the gluon and decoupling of the gluon polarization effects. 
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Figure 1 : The running coupling of the do- 
main wall fermion. Coulomb gauge gluon- 
ghost coupling of m u = 0.01/a(square), 
0.02/a(diamond), and quark-gluon cou- 
pling of m u = 0.01/a(large disks). Small 
disks are the a s , gi derived from the 
spin structure function of the JLab 
group|33]and the solid curve is their fit. 
(color online) 
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Figure 2: The beta function calculated 
from the logarithmic derivative of the run- 
ning coupling parametrized by the JLab 
group. 



3. Schrodinger functional and the critical flavor number 



Recently Appelquist [13|, [14j claimed by performing lattice simulation of 



running coupling in the Schrodinger functional method [171 . [18[ that there is 
a critical number of flavors Nf c below which both chiral symmetry breaking 
and confinement set in. He assigned 8 < Nj < 12 and in the case of Nf = 8 
the running coupling monotonically increase as decreases and in the case 
of 16 4 lattice and f3 ~ 4.65, g 2 (L)/An ~ 20/4tt ~ 1.7, and that there is no 
sign of infrared fixed point. However, the Schrodinger functional is derived 
through the field that satisfies a boundary condition at x = and x = T 

B = 0, B k = [x°C' k + (T-x°)C k ]/T (2) 

where in the case of SU(3) where 

/ <j> kl \ / <j> hV \ 

C k = - fc2 andC^ = - fe2 , (3) 

V 0*3 / V ) 

in which conditions ^ = and ^^^ki' = for all k is adopted. 

The functions <fi and <$' are parametrized by rj and v and the running 
coupling is derived from 

^=^--r (4) 
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Q 

The Schrodinger functional method a s (q) is parametrized as 



logg/A 

and its applicability in the region of q below and around A ~ 213 ± 40MeV 
is not clear, since the vector field could be defined from fermionic field and 
the structure of the manifold above and below A could be different. 

Brodsky et al, argue that in the DSE q < A or r > lfm, and in 

ADS/QCD ( = y/b±x(l — x) > A 

max ~ lfm, confinement is essential and 
the quark-gluon should not be treated as free. 

The QCD coupling constant is related to the bare coupling as 

9 = Z//*Z^g (5) 

and in one loop perturbation theory, the gluon self energy yields 



a 2 13 

1 + 0*200 logA//. (6) 

lD7r 3 



and the vertex renormalization factor is 

where £ = 1 in Coulomb gauge and 

t 2 (V)5 ab = C acd C bcd (8) 

which is n5 ab for SU(n). Our conjecture is that, in the infrared region, 
the vector potential is produced as Pliicker coordinate of fermions and that 
triality symmetry of the fermion effectively triplicate t2(V). 



4. Anomaly matching condition and number of flavors 

Infrared QCD is characterized by spontaneous chiral symmetry breaking. 
The quark field of QCD at low energy is expressed as[l| 

l/> = -(l+J 5 )ll>L+^(l-J 5 )ll>R (9) 

where ipL transforms as SU(3) color under G c , SU(N) L flavor under Gp and 
SL(2, c) Lorentz spinor and transforms as SU(3) color, SU(2) R flavor 
and SL(2, c) Lorentz spinor. 
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Color gauge fields bind these quarks into baryons which must be color 
singlets. In the model that gauge group G c coupled to chiral fermion, the 
representation of the group must be such that the anomaly is canceled. In 
the space of symmetry G c x Gp, t'Hooft considered 

G F = SU{ ni ) L g> SU{n 2 ) R ® SU{n 3 ) L <g> SU(n 4 ) R g> U(lf (10) 

where n 12 refer to the triplets and n 3 4 to sextets. He proposed anomaly 
matching relation in which the symmetric anomaly coefficient d abc contributed 
from the massless quarks should be equal to the coefficient D abc contributed 
from the massless color singlet composite fermions jlj]. 

The approach up to fermionic degrees of freedom Nf = 5 was not suc- 
cessful, but an algebraic research was extensively done in jij , and it is shown 
that the matching of SU(Nf) 2 L x U(1) B and SU(Nf) 2 R x U(1) B anomalies, 
where U B (1) is the vector-like Baryon number, and matching of SU(Nf) 3 L 
and SU(Nf) R anomalies can be done in Nf ± 6 dimensional representation 
0| . The condition on SU(Nf) 3 L and SU(Nf) R anomalies is 

3 = -(N f ± 6) + i(JV> ± 3)(JV> ± 6) + l -N f {Nf ± 1) - N f (N f ± 4) (11) 

and that for SU(Nff L U B (l) and SU(N f ) 2 R U B (l) anomalies is 

1 = -(Nf±6)^^ + ^(Nf±2)(Nf±3) + ^Nf(Nf±l)-Nf(N f ±2). (12) 

where U B (1) is the baryon number charge. It means that QCD with 3 flavor 
and 3 color system cannot satisfy the matching condition. 

Sannino [4] found another solution that satisfies the anomaly matching 
condition in 2Nf ± 15 dimensional system, which requires Nf > 8. He found 
the condition for a SUlNf) 3 ^ anomalies as 

3 = -{2N f ± 15) + \{N f ± 3)(JV) ± 6) + l -N f {N f ± 1) - N f (N f ± 4) (13) 

and that for SU(Nf)j j Uv(l) anomalies as 

1 = -{2N f ± 15) ™ f ±* 5 +\(N f ± 2) (N f ± 3) + ±N f (N f ± 1) - N f (N f ± 2) 

(14) 

where Uy(l) is the vector charge. 
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The critical number of flavors for opning the conformal window [9] and the 
critical number of flavors for presence of infrared fixed point is expected to 
be close to each other. A study of infrared fixed point was done on the lattice 
by measuring the momentum dependence of the QCD running coupling. 

In the perturbative /3 function and DSE approach, infrared fixed point 
am is a decreasing function of number of flavors Nf. When am decreases 
below the minimal value of a cr for which the fermion acquires a mass, the 
Nf is called Nf c 16l. Il9j]. Near this Nf, the effective coupling varies slowly as 
the momentum changes. Lattice simulation of N c = 3 [5] suggests Nf ~ 7, 
while more recent lattice simulation [Lj, [14j, |l5| and DSE[16| suggest that 
8 < N f c < 12. 



5. Octonion and Anomaly matching condition 

In the vector like theory that SU(N f ) 2 L x U{l) v and SU(N f ) 2 R x U(l)v 
anomalies and SU(Nf)^ and SU(Nf) 3 R anomalies can be matched in 2Nf ± 
5N C dimensional representation The model of j3] may be regarded as 
magnetic dual of the gauge theory. The factor 5 in this QCD dual may 
be related to the 1/N C corrections to the one body matrix elements in the 
Skyrme model [4l| which was necessary to reproduce the magnetic transition 
of a nucleon to A. These condition Nf — 6 = 3 and 2Nf — 15 = 3 are both 
satisfied by Nf = 9, but the standard model of QCD predicts Nf = 3. The 
quadratic form of quaternion is expressed by octonion, and the automorphism 
in the space of octonion is G 2 group which is exceptional Lie group with 14 
dimensional representation. In G2, there is a specific automorphism, which 
is called triality. In the following subsection, I study the structure of the G2 
group. 

5.1. Representation theory of G2 group 

In this subsection, I study the quaternion and the representation theory 
of (72 0. 

A quaternion geH can be defined by using the basis {1, i, j, k} = {I, io~\, io~2, «cr 3 } 
as q = wl + ixo~i + iyo~2 + iza 3 , where <Tj are Pauli matrices and / = l 2 is the 
unit matrix. As elements of a complex number is expressed by a product of 
real numbers C = 1Z + ilZ: 

(xi, y l ){x 2) y 2 ) = (x 1 x 2 - ym, xm + yix 2 ) (15) 
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a quaternion is given from bijection of two complex numbers 

Ci 



*- 


2 ■ 


*-( 


Z 2 


w 2 \ 






v -w 2 


Z 2 ) 



and the elements of quaternions are expressed by a product of complex num- 
bers as the first low of the product of £i and ( 2 shows: 

{z 1 ,w 1 )(z2,w 2 ) = (z 1 z 2 - Wxw^ z x w 2 + Wiz 2 ) (16) 



Thus "H = C+jC |43j. The quaternion H and a new imaginary unit 1 that 
anti-commute with i,j,k compose an octonion O = H + YH. It is spanned 
by 

{1, i, j, k, 1, il, jl, kl} = {1, ei, e 2 , e 3 , e 4 , e 5 , e 6 , e 7 } (17) 

with bases of Clifford algebra O = R + R 7 . 

Elements of an octonion is expressed as a product of quaternions 

(pi, qi) (P2, 92) = (P1P2 - Mi, Q2P1 + qm) (18) 



The algebra O is associated with an involution [44 

9^9= (I ) ^ ( -°i) ^ 

where 7 1 is the 7 dimensional diagonal matrix. In the subspace R 7 , automor- 
phisms of 50(7) are of the form U — > SU S^ 1 where SeSO(7). Automorphism 
of Spin(7) are of the form u — > sus -1 . 

In the space R 8 , automorphism of Spin(8) is expressed as injection 

H'-mS):.^ , v ° {a) ) (20) 

whose image is defined as Y. For an element y G Y and the unit element 
g e O, involution g of g is defined as 



gye = gye. (21) 

where ye is the complex conjugation of ye. When the element g satisfies 
ge = e and g also satisfies ge = e, g is called the companion element of g. 
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In the case of SO (8), with g and the companion of g%, an element of 

7 \ 

) and by the definition of g 2 of 
u 9i J 

goVSlT 1 * = faye (22) 

for all y G Y, one can construct a triple ( y g ,gi,g 2 ). 
The triality automorphism is 

e:Spin(8)^Spin(8);( 9 ° ° ) ^ ( q M (23) 



Triality transformation rotates 24 dimensional bases defined by Cartan|45 



{6, 6, 6, 6, 6}, {62, 6l, 63, Cl4, ^24, £34}, {£l23, 624, 6l4, 634, 6234} 
r 1 2 3 4i r 1' 2' 3' 4'i 

J /■>' T* T T* > J ">"■ r f* rf T* > 

The trilinear form in these bases is 

J= = <t> T CXilj = x l (i l2 £- ilA - £31624 - 614623 + 62346) 

+ ^ 2 (63624 — 62634 — 64623 + 62346) 

+ a; 3 (61634 — 63614 — 64623 + 62346) 

+ x 4 ("64634 — 64614 — 64624 + 62346) 

+ x 1 ("6634 + 636 — 646 + 646) 

+ a; 2 ' (-6614 + 616 - 646 + 646) 

+ a; 3 ' (-6624 + 626 - 646 + 646) 

+ x 4 ' (£0623 - 636 - 616 - 626) (24) 

There are three semi-spinors which have a quadratic form which is invariant 
with respect to the group of rotation 

$=*0C0 = 66234-6364-6164-6264 (25) 
tf=*^Cty = 6634-6634-6624-6623 (26) 

and the vector 

F = x 1 x 1 ' + x 2 x 2 ' + x 3 x 3 ' + x 4 x 4 ' (27) 

The triality transformation that makes ge = e for e Ob >j Ob « 00 and 

x x \x 2 \x 3 ' is defined as G 23 , ge = e for e = 6,6,6 and 624,614,634 is 
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defined as Gi 2 , ge = e for e = £12, £31, £23 and £14, £24, £34 is defined as Gi 3 . 
In additioin to G 2 3,Gi 2 and Gi 3 , there are automorphism t (G\ 3 Gi2) = Gi 32 
and i {Gi 2 G 13 ) = Gi23. 

When one defines for each i = 0, 1, 2 

Hi = {(90,91,92) e Spin(8) : &e = e}, 
the triality transformation is expressed as 9(Hq) = H, 0(Hi) = H 2 and 



9{H 2 ) = H , and G 2 = H x n H 2 = H 2 n H = H n #i|44|. 

In the space of octonions, I choose ^ as an element orthogonal to 1, and 
£2 as an element orthogonal to 1 and £1, and £4 as an element orthogonal to 
1, £1,^2 and £i£2. For orthogonal elements £i,£2,£4 chosen as £&, there is an 
automorphism £1 — » £14, £2 — > £24, £4 — > £o- The choice of £14 is a choice of a 
point on a unit sphere in the 7 dimensional space, £24 is a choice of a point 
on a unit sphere in the 6 dimensional space, and £ is a choice in the space 
orthogonal to 1, £ i4 , £24, £u£ 2 4 i-6. in the 4 dimensional space. These form 6, 
5 and 3 dimensional manifolds. The exceptional Lie group representation G 2 
is the automorphism in this 14 dimensional manifold. 

I assign SU(3) fundamental color representation of quarks on the 3 di- 
mensional manifolds and that of the color representation of diquarks on the 6 
dimensional manifolds and 4 dimensional spinors and 1 dimensional space be- 
tween spins on the domain walls to the 5 dimensional manifold. The triality 
transformation mixes the 3 and 6 dimensional spaces. 

Irreducible representations of g 2 algebra with highest weight auj\ + bu 2 is 
written as r a b , and I study the standard representation V = T 1Q . 

5.2. Anomaly matching condition 

A problem is whether the anomaly matching of the color triplet and 
color sextet sectors in a baryon can be explicitly realized. The fundamental 
representation of the g 2 algebra has 14 vectors 

92 = {Hi, H 2 , Xi, Yi, X 2 , Y 2 , X 3 , Y 3 , X 4 , Y 4 , X 5 , Y 5 , X 6 , Y 6 } (28) 

One can define a subgroup 

g = {H 5 , H 2 , X 5 , Y 5 , X 2 , Y 2 , X 6 , Y 6 }, (29) 

where H 5 = Hi + H 2 , which is isomorphic to sl^(C) (42|. 
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The rest of the Lie algebra consists of 

W = {X,,Y 1 ,Y 3 } and W* = {Y 4 , X u X 3 } (30) 

When a representation has the highest weight uj\ = 2«i + it is seven 
dimensional, we call it V . The wedge product A 2 V is 

A 2 v ~ r 0j i + V (31) 

and A 3 V is 

A 3 V ~ r 2i0 + V + C (32) 

where C is the trivial representation. 

The action of gi on the standard representation preserves a skew-symmetric 
trilinear form u> on V, which consists of 

oj = w 3 A u A v 3 + v 4 A ii A w 4 + Wi A w A v i + 2^i A f 3 A w 4 + 2w\ A w 3 A v 4 (33) 

where with a use of root vectors Y\ and F 2 , ^3 = ^1(^4), ^1 = — ^2(^3), w = 
Yi(t>i),u>i = Y 1 (-u)/2,w 3 = 12(101) and u4 = — Fx(w3). It means that the 
five- dimensional space u is stable in A 3 V, and that the anomaly matching 
relation of the massless fermion state and the bound state of three massless 
fermion which is given by A 3 V can be considered in the same five dimensional 
space. 

sw(3) is a subalgebra of s/ 3 (C), but it is a real representation. The quater- 
nion real condition imposed on the quark propagator has the effect of con- 
straining sls(C) to its real representation su(3). Anomaly matching relation 
on the massless fermion and the bound state of three massless fermions with 
flavor Nf = 3 and color N c = 3 could be interpreted as the Nj = 9 without 
triality, since the same 5 dimensional space appears 3 times through triality 
transformation. The topology around the infrared fixed point would be more 
complicated than the U(l) 3 . 



6. Gauge fixing of the domain wall fermion 

We first perform a gauge fixing in Landau gauge using the over relax- 
ation method and then fix to the Coulomb gauge without touching the A 4 



component [471] . I then measure the quark propagator and fix the gauge in 



the space of the 5th dimension between the domain walls. 
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In the analysis of instantons in quaternion bases, Corrigan and Goddard 
46[ defined the transition function g(u>, it) where u and 7r are complex two- 



spinor which satisfy 

g(\u, \n) = g(u,7i), detg = l. 
©When u = XTi where x = x —ix-cr, g(u, 7r) is a quaternion can be expressed 



as 



g{x7c,ir) = h(x,()k( x >0 1 
where ( = — , h(x,() is regular in \(\ > 1 — e and k(x,() is regular in 

ICI < 1 + c n 

In 46 , the Ansatz 



C 1 P \ ( e M 



• ,/ " 1 o e ^ y v c 1 y v ° ' " 

e^C 1 



(34) 



was proposed as the transformation matrix. In our 5-dimesional domain wall 
fermion case, 7 = /j — v and //, z/ contain the phase in the 5th direction irj. 

2/i = iLOij'Ki — ir\ = (xi + ix2)( + ixo — £3 — ir] (35) 

2z/ = iu)\/it\ + ir] = (x± — ix 2 )( + ^0 + £3 + iv (36) 
The quaternion reality condition of the transformation matrix (7(7, 



gives 



A) 



(37) 



where / = /(7,-j). 

The function / and / taken in |46( is 



doe 7 — e 7 -3—^ — e 7 — ane 7 

yr _ a.L s -i j _ dj^-i u 

In general c and cl s _i are polynomials of (, and they satisfy c^-iC 1 
CoC -1 - I define 

V' = c„ 1 C" 1 + c 1 C 1 + 5 
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where c_i = c^ s _i and C\ = c and 5 is a constant, which is defined later. 
With the Ansatz ip = c C _1 + c L ^C 1 + 5, I find C 1 = J and the e 7 

\l c L 8 -l 

and the 5 are calculated from the simultaneus equation (j6]): 

f a^-i/ + K-iC'e^ = koC^e" 7 + d f 
\ c Ls .J + d L .-xC- 1 e- rl = d C l e\ 

The equation has two sets of solutions, which can be obtained numerically. I 
assign ( a ° J° j and ( aLs ~ x \ Ls l J from the color diagonal com- 

V Co d J L/R V C Ls-l 0-L.-1 J L/R 

ponents of the lattice data. 

From the numerical practice, I observe that the deviation AL/R becomes 
small when \5\ is large, i.e. when ip ~ 5 ~ j is large. It means that 

e 2 7 1/8-1 is a good approximation in the case of small deviation, and 
d 

I select the gauge such that the corresponding \5\ is large. 

7. Form factor of the proton 

With the prescription on the gauge of fermion on the domain walls, I 
calculate the baryon charge form factor of a proton using the SU(6) spin- 
flavor wave function. The nucleon three point function is 

(G^(t 2 ,t i;P + |,p-|;r)) 

= ^ e -(P+f)-*v(P-f)-*a x r ^<n|T|x a (x a) t 2 )^(o,o)^(a! 1> -t a )]n> 

x 

where 

x (x) = e abc [u Ta (x)C l5 d\x)]u c (x) 

In these expressions abc specify colors, a and (3 specify a spin dependent A 
term or a spin independent B term of each three quark propagators. 

In the limit of t\ — > —00,^2 — > +00, the fourier transform reduces to the 
vertex sandwiched by the propagators. 
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7.1. The case of q = 

I define the normalization by the ratio of the three point function G NVflN 
and the two-point function defined as 

(G NN {t2,t u p;T)) = X)e^ aa e^r^<n|Tx a (a! a ,t a )^(a!i,ti)|n>. 

x 

I take the limit of t\ — > — oo, £2 - ^ +00, here also. I define the propagator 
S(a,h a ,a,h' a \p) = S A (a,h a ,a,h' a \p) + S B (a,h a ,a,h' a \p), (38) 

where 

S A (a,ha,a,h'M = Ap ~^ P MB \KK 

B 

S B (a, h a , a, h' a \p) = Ap2 + MB \kk 

In the coupling of A M to (u a (x)C'~f 5 d b (x))u c (x) and (u a ' \x)C"y 5 d b ' (x))u c ' (x), 
I evaluate the real part of 

1. S a (c, h c , c, h' c \p)u\(p)A (t u c h » c (p) x lhS p {c, K\, c, /i c |p) f 7 5 

2. S Q (a, h a , a, h' a \p)u a K (jp)K»u a h » a (p) x 7s ^(a, h\, a, fr a |p)t 7 5 

3. ^(6,-^,6,-^1 -p)^u b K {-p)k„u\ % {-p) x ^(6,-/1%, 6, -/i 6 | - 

I used the fact that the helicity of C^d{x) becomes -1 times the original. 
The expectation values of the charge form factor G^+ N is 

J 2 [ Gf/(p,0) + e>,0) + G™( P ,0) + G™(p,q) 

,fVM , G NVN b+ ( P ,o) 

+ G^V(p,0) + G^ 6+ (P,0) J 

At zero momentum the left-handed quark contribution dominates when the 
gauge fixing parity on the domain wall fermion is applied. When the gauge 
fixing parity is not applied, the left-handed quark and the right handed quark 
give almost the same contribution. 

I show in Fig. [21 the form factor of the proton of the domain wall fermion 
of 16 3 x 32 x 16 lattice of ud— quark mass 0.01/a (148 samples). The momen- 
tum p is chosen to be (pi,p 2 ,p 3) P4) = (0, 0, 0, 0), (1, 0, 0, 0), (1, 1, 1, 0), (4, 0, 0, 0), 
(2, 2, 2,0), (3, 3, 3,0), (4, 4, 4,0). 

The dotted line is the dipole fit with M 2 = (0.71GeV) 2 . Since the con- 
tribution of left-handed and right-handed are almost the same for p ^ the 
fitted line is normalized to be 0.5 at p = 0. 
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p(GeV) 



Figure 3: The form factor of a proton using the DWF m u = 0.01 /a gauge fixed such 
that the fermion on the left wall and the right wall arc correlated by self-dual gauge field, 
and the dipole fit with M 2 = (0.71GeV) 2 . At zero momentum left-handed(LH) fermion 
dominates the form factor, and at finite momentum, LH and RH contribute almost equally, 
(color online) 



7.2. The case of q ^ 

In 31| . the local vector current Ward identity on quark state is written 

as 

q^Zythy,^ + \ V ~ §)] = -*(S-\p + |) - S-\p - |) (39) 
q q 

where A V L^(p + -,p — — ) is the amputated vertex with momentum transfer 

q and Z v l is the renormalization factor. When one operates d/dq p on this 
expression, one obtains 

Z vL (A L Vp (p) + q^K P (P + %V ~ §)U) = -^^(P) ( 4 °) 

q q 

In our lattice simulation, I take the difference of vertex of p H — and p 

F 2 F 2 

in which q is parallel to p, I consider the difference near the cylinder cut 

----- Q 
momenta p = (p,p,p,2p), (p = 1,2,3,4), i.e. the difference of p — - = 

(1, 2, 2, 4), (2, 3, 3, 6) and P+tj = (3> 2, 2, 4), (4, 3, 3, 6), respectively in the case 

of p = 1, and p - | = (2, 1, 2, 4), (3, 2, 3, 6) and p + | = (2, 3, 2, 4), (3, 4, 3, 6), 
respectively in the case of p = 2 . 

I calculate for each three body vertex G^ VN (p, q) in S^{p ± -) as 

(u a (x) C^d\x))u c (x) | S^A^ | (w a (x)C7 5 d b (x) ) w c (x) ) 
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Cuuaa Cuua C-uubc Cpp C-pu C-pua C-up C-upa 



4/18 1/18 -2/18 1/18 1/18 -2/18 1/18 -2/18 
1/18 4/18 -2/18 1/18 1/18 -2/18 1/18 -2/18 



= S(ca,ca'\p + ^)[u c a'A fl u c a"] l5 S(ccr",c ( T\p+ |)t (41) 

Although gluon can change color of the quark, I select color diagonal com- 
ponent and spin off-diagonal components. 

The nucleon helicity a is chosen to be + and but at the vertex A, the 
helicity of the quark may change through the propagators'. 

The two body part G NN (p ± ^) is evaluated by making a scalar product 

of 

/ T$ a aa „ \ / Abs[S(aa', aa» \p ± §) £ + S(aa', aa" \p ± §) fl ]/2 \ 
T= T <U" } = [ Abs[S(ba' b ,ba"\p±l) L + S(ba' b ,b<j"\p±l) R }/2 
\ T c a , (j c " / V Abs[S(ca>, ca» \p ± f ) L + S(ca' c , ca» \p ± § ) R }/2 J 

(42) 

where a' a" etc are ++, H — , — h and , and 

/ S bc (a b a c ,a' b a c ») \ j S(ba b ,ba' b \p + f) L x S(ca c ,ca c "\p + f)^ \ 
S= S ca (a c a a ,a' c a a ») = S( C(7c , c<|p + §)l x S(aa a , aa a " |p + f)^ 
V S a6 ((7 a( 7 6 , <o- 6 ") / V 5(aa a , aa> + \) L x S(K, W \p + l) R I 

(43) 



N=(<Zoto*\T.S\o*of<t) 

where a a ai,cr c are + H — , H h or — h + and according to whether the flavor 

of the quarks with helicity ++ coincides or not, the normalization of the 

• 2 1 
three quark state is — = or — . 

VT8 VT8 

The three-body vertex matrix elements G NVN (p ± ^) is calculated as 

follows. I define the vertex matrix element on the left-handed quark 'a' and 
spectator quarks 'b' and 'c' as QxyLXz where xy specify whether the helicity 
of spectators in the initial state and the final state respectively are parallel(p) 
or unparallel(u), z defines the type of helicity flip of the spectators and X 
defines the type of helicity flip of the vertex. The corresponding two body 
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matrix element is defined as NxyLXz. In the calculation of left-handed 
quark contribution in A a / aa » , I measure 



Q = tfoto*\Sh& ■ S\a*a*a L c ) 



(44) 



where 



S(aa a , aa'Jp + Z) L A <<Ta »S(acr a , aa a " \p + 
S(ba b , ba' b \p + Z) L A a > ab »S(ba b , ba b " \p + 



2>R 
2>R 



S(ca c , ca' c \p + |) L A a / CTc "S'(co- c , «r c " \p + |)Jj 

(45) 

as a sum of the matrix elements given in the Appendix. 
I define 



AppLl ■■ 
AppL3 - 

ApuLl ■■ 
ApuL3 ■ 
AupLl ■■ 
AupL3 ■ 
AuuLl 
AuuLI 
AuuL3 
AuuLA 



Re[QppLl}/Re[NppLl], AppL2 = Re[QppL2}/ Re[NppL2], 
Re [QppL3] I Re \NppL3\ , AppLA = Re [QppLA] /Re [NppLA] 

ApuL2 = 
ApuLA = 
AupL2 = 



Re 


Cp U QpuLl + 


Cp Ua QpuLla] 


Re\ 
Re 


CpuNpuLl + 
c pu QpuL3 + 


Cp Ua NpuLla\ 
Cp Ua QpuL3a} 


Re\ 
Re 


c pu NpuL3 + 
c up QupLl + 


Cp Ua NpuL3a] 
c upa QupLla] 


Re\ 
Re 


c up NupLl + 
c up QupL3 + 


c upa NupLla] 
c U p a QupL3a] 


Re\ 


c up NupL3 + 


c upa NupL3a] 



AupLA 



Re[c pu QpuL2 + 


Cp Ua QpuL2a] 


Re[c pu NpuL2 + 
Re[c pu QpuLA + 


Cp Ua NpuL2a] 
c pua QpuLAa] 


Re[c pu NpuLA + 
Re[c up QupL2 + 


Cp Ua NpuLAa] 
c upa QupL2a\ 


Re[c up NupL2 + 
Re[c up QupLA + 


c upa NupL2a] 
c upa QupLAa\ 


Re[c up NupLA + 


c upa NupLAa] 



Re[c uua QuuLla + c uuaa QuuLlaa + c uubc (QuuLlb + QuuLlc)] 
Re[c uua NuuLla + c uuaa 
Re[c uua QuuL2a + c 



NuuLlaa + c uubc (NuuLlb + NuuLlc)] 
QuuL2aa + c uubc (QuuL2b + QuuL2c)} 



(46) 



Re[c uua N uuL2a + c uuaa NuuL2aa + c uubc (NuuL2b + NuuL2c)] ' 
Re[c uua QuuL3a + c uuaa QuuL3aa + c uubc {QuuL3b + QwuL3c)] 



Re[c uua NuuL3a + c u 
-Re[c uua QwwL4a + 



1 NuuL3aa + c uubc (NuuL3b + A/"miL3c)] ' 
a QuuLAaa + c uubc (QuuLAb + (3uwl<4c)] 



(47) 



Re[c uua NuuLAa + c uuaa NuuLAaa + c uubc (NuuLAb + Af-uwL4c)] 
and diagonalize the left-handed quark contribution 

/ AppLl + ApuLl + AupLl + AuuLl AppLA + ApuLA + AupLA + AuuLA 
I AppL3 + ApwL3 + AwpL3 + AtmL3 AppL2 + ApwL2 + AupL2 + AuuL2 
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Figure 4: The charge form factor of a proton using the DWF m u = 0.01/a gauge fixed 
such that the fermion on the left wall and the right wall arc correlated by self-dual gauge 
held, and the dipole fit with M 2 = (0.71GeV) 2 . (color online) 




take eigenvalues and multiply the normalization factor 18 2 . I calculate also 
the corresponding matrix elements of right-handed fermions, and add the 
two. 

Number of samples is 52 and the error estimation was done in by using 
Bootstrap resampling command of Mathematica with 5000 iterations. (I 
calculate standard deviation of left-handed quark contribution divided by 
the mean of the left-handed quark contribution and the corresponding value 
of the right-handed quark contribution and multiplied the mean to the sum of 
the left-handed quark contribution and the right-handed quark contribution.) 

In contrast to q = case, the fluctuation is relatively large, as shown in 
the FigJU The contribution of Ay^blue (dark)) and Ay 2 (green) are consistent 
with a dipole fit of M 2 = (0.71GeV) 2 within errors. The standard deviation 
of Ay 2 is larger than that of Ay 2 . 

The nucleon electromagnetic form factor of the domain wall fermion was 



calculated recently by LHPC collaboration [50[. The electric form factor 
Ge{q 2 ) in their simulation has a small pion mass dependence, but the mo- 
mentum dependence is flatter and the absolute value is larger than the phe- 
nomenological dipole fit. I guess the main difference from this calculation 
comes from the topology of the manifold. In our simulation, the helicity 
hih 2 of \hih2 an d T hlfl2 are the same, but when the boundary condition 
of the three quarks of the vertex QxyLXz and that of the normalization 
NxyLXz are fixed to be the same, hih 2 of the two are not necessarily the 
same and in this case the momentum dependence becomes flatter. 
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8. Discussion and Conclusion 



I showed that the discrepancy on the critical flavor number of fermions Nj 
for the appearance of the infrared fixed point derived from the Schrodinger 
functional method and the Bjorken sum rule and GDH sum rule and our 
lattice simulation could be resolved by assuming the topological structure 
of the infrared fixed point is not that of U(l) 3 but that of G2 with triality 
automorphism. 

The Lie algebra gi has the 8 dimensional subspace go which is isomorphic 
to sl 3 (C) and the 3 dimensional subspace W and W*. si 3 (C) is not iso- 
morphic to su(3), but by imposing quaternion real condition, the subspace 
isomorphic to su(3) could be selected. 

The lattice simulation of QCD running coupling from quark-gluon cou- 
pling in Coulomb gauge suggests infrared fixed point of a s (0) ~ n which is 
predicted in the BLM renormalization scheme. The strong interaction sum 
rules obtained in the infinite momentum frame with holographic variable z in 
addition to the 3+1 dimensional space time is 5 dimensional and similar to 
ours. I expect the similarity originates from the similar topological structure 
of the manifold, i.e. S 3 x R of the quaternion. 

Effects of instanton on S 3 x R topology in Coulomb gauge was studied in 
the quaternion bases [ill]. I observed that the self-dual gauge field (instanton) 
plays important roles in the stability of the infrared QCD. 

Using the domain wall fermion in quaternion bases, and taking into ac- 
count that the quadratic term of a quaternion makes an octonion, I remarked 
that there are 3,6 and 5 dimensional stable manifolds and the exceptional 
Lie algebra is the automorphism in this space. The gi has triality au- 
tomorphism and the topology of the infrared fixed point is not trivial. I 



extended [48[ and [49J to a larger number of gauge configurations, and cal- 
culated the form factor at q = of the proton choosing a specific gauge that 
the fermion on the left domain wall and the right domain wall are correlated 
by the instanton-like self-dual gauge fields. I observed also via Ward identity 
through the difference of the inverse propagator at momentum p+q/2 and 
at p-q/2, charge form factor of the proton at q ~ 0.4GeV and 0.9GeV agrees 
with the phenomenological dipole fit. Our simulation differ from the standard 
lattice simulation due to different boundary condition on three quarks. 

The triality automorphism gives another possible solution to the U(l) 
problem, i.e. the u^— , d^— and s^— quarks need not have the same phase e~ 10 
and the ur— , d R — and sr— quarks the same phase e ld due to extra degrees 
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of freedom. It could affect or change the interpretation of the width of 
i] — > 7r + 7r°7r~ decay, which is about factor three larger than the standard 



current algebra result [53 



I used information of the color diagonal matrix elements only of quark 



propagators, and ignored the vector product part of quaternions 43] . I need 



to calculate other form factors of a proton and check the correlation of the 
three quarks in a proton. It is necessary to extend the simulation to a larger 
lattice and calculate form factors at more momentum points and examine 
the triality automorphisms in more detail. 



I would like to thank Dr. Francesco Sannino for valuable information on 
the ano maly matching, Professor Stan Brodsky for attracting attention to 
the ref. [l9|] and helpful discussion, and Dr. Alex Deur for sending me the 
experimental data. I acknowledge helpful suggestion of paying attention on 
recent DSE approach [7J and on the Pinch Technique [a], from one of the 
referees of a journal. 

The numerical simulation was performed on Hitachi-SRllOOO at High En- 
ergy Accelerator Research Organization(KEK) under a support of its Large 
Scale Simulation Program (No.07-04 and No.08-01), and on NEC-SX8 at 
Yukawa institute of theoretical physics of Kyoto University. 
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Appendix .1. 7 matrices and triality 
I use the 7 matrices, as follows: 




75 



/ 1 \ 

10 

0-10 

\ -I ) 



74 



/ 1 \ 

1 

10 

\ 10 0/ 

(•1) 



71 = 



/ 











1 \ 




/ 











-1 \ 




/ 








i 





\ 










i 















1 


















—i 









—i 








72 = 







1 








73 = 




—i 













V 


—% 












V 


-1 








) 




V 





% 








/ 



(•2) 







( 


1 





o\ 




/ 





1 








\ 








-1 















-1 















747275 = 













1 


7472 = 













-1 










V 





-1 


0/ 




V 








1 





/ 






/ 





—% 





N 






/ 


-% 








o\ 






—i 




















i 










747i75 = 













—i 


747375 = 












— i 









V 








—i 


) 






V 











i 


I 



(.3) 



(•4) 



22 



Appendix .2. The matrix elements of the triality transformation of the 
algebra. 

In G 23 , I replaced £ 3 -»■ £ 34 to £ 3 -»■ -f 34 , and inGu, £ 4 -> -£1234 to 
£4 — >■ —^123 and x 4 — £ to s 4 — y £ n written in [45] (probably typos). 
G13G12 does not agree with Gr\ss\ of [45]. I think £ 4 — > —x 4 there should 
be replaced by £4 — > x 4 and I define t (G 12 G r i 3 ) = G123 and t (G 13 G 12) = 
G132, which correspond to G(i23) and G(i 3 2) respectively. Under this triality 
automorphism, the set {F, <3>, is invariant a signe pres. (In the tables 1 
indicates -1.) 
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Appendix .3. The matrix elements of the vertex QxyLXz and the normal- 
ization NxyLXz in the form factor calculation 

In the calculation of the G^ VN with A M=1 , I used the martix elements 
of the vertex QxyLXz and the normalization NxyLXz as follows and the 
corresponding matrix elements of L and R interchanged, i.e. QxyRXz and 
NxyRXz. The vertices of G^ VN with A M=2 is also calculated similarly. 
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NuuLAc = l{T-+ • S-+-+)r, 

QpuLl = L { - + + | ( 5_+A ++ ^ + • S+++4 + +-) R , 
NpuLl = L (T ++ ■ S +++ -) R , 

QpuLla = L ( - + + \S-+A ++ Sl+ ■ + -+) R , 

NpuLla = L (T++ ■ S ++ _ + )r, 

QpuLl = L ( - + + |S__A__Sl + ■ + +-) R , 
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QpuL2a = L { - + + |S__A__SL + ■ S ++ - + \ + - 
NpuLla = L (T__- S ++ - + } R , 

QpuL3 = L ( - + + |S_+A+_Sl + ■ S +++ _\ + +-) R , 
NpuL3 = l(T+- ■ S +++ _)r, 
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QpuLA = L ( - + + IS 1 — A_+5'| + • + +-) R , 
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NpuLAa = L (T-+ ■ S ++ _ + ) R , 
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NpuLl = L (T ++ ■ S + - ++ ) R , 

QupLla = L ( + - + \S ++ A ++ Sl_ ■ S_ +++ \ - ++) R , 
NupLla = L (T++ ■ S_ +++ ) R , 

QupLl = L ( + + - \S + _A__Sl_ ■ S + _ ++ \ - ++) R , 
NupL2 = L (T— ■ S + - ++ ) R , 

QupLla = L ( + - + \S+-A—Sl_ ■ - ++) R , 

NupLla = l(T__- S- +++ ) r , 

QupL3 = L ( + + - \S ++ A + _Sl_ ■ S + _ + +\ - ++) R , 
NupL3 = L (T + _ • S+- ++ ) R , 
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NupL4 = L (T_ + -S + _ ++ ) R , 

QupLAa = L { + - + \S + _A_ + Sl_ ■ S_ +++ \ - ++) R , 
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QppLl = L ( - + + \S- + A ++ Sl_ ■ S ++++ \ - ++) R , 
NppLl = L (T ++ -S ++++ ) Rl 

QppLl = L ( - + + |S__A__Sl_ • S ++++ \ - ++) R , 
NppL2 = L (T__-S ++++ ) R , 



Q PP L3 = L ( - + + |5 , _+A+_5l_ ■ - ++) R , 
NppL3 = • S ++++ ) R , 

QppLA = L { - + + |S__A_+S|_ ■ - ++} R , 
NppLA = L {T_ + -S ++++ ) R . 
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